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Abstract 

Some results for commuting local integrals of motion for free bosonic theories are presented, which 
can be applied in construction of integrable boundary states. Usable program code for calculation 
of these integrals of motion is also introduced. 

1 Introduction 

In the set of theories in two dimensions, conformal ones are fixed points |^j with respect to Renor- 
malization Group (RG) transformations |19j . Perturbing of such a theory with certain relevant 
operator yields integrable flow, i.e. a trajectory of "deformed" theories, which leads to another 
fixed point, which is usually called "infrared fixed point". The intermediate theories on such a 
trajectory are integrable, meaning that a countable subset of the integrals of motion (IM) from the 
original theory (which is sometimes called "ultraviolet fixed point" ) survive the relevant perturba- 
tion. In the last several decades there has been a substantial progress in the understanding of the 
relation between Conformal Field Theory (CFT) and Integrable Quantum Field Theory (IQFT) 
in two dimensions [T] |5] [3] |£] [TT] [T71]. 

A conformal field theory is characterized by its chiral algebra A. The algebra A can be either the 
Virasoro algebra or an extended algebra which includes the Virasoro algebra as a subalgebra. The 
Virasoro algebra is generated by the holomorphic component T of the energy-momentum tensor. 
The antiholomorphic component of the energy-momentum tensor T correspondingly generates 
another (antiholomorphic) copy of the Virasoro algebra. From now on we will restrict our attention 
to the holomorphic part of the theory, assuming that the theory possesses reflection symmetry and 
therefore all the reasoning about the holomorphic sector pertains to the antiholomorphic one as 
well. 

One of the fundamental characteristics of an integrable field theory is the presence of a countable 
set of commuting integrals of motion {I s }, 

[i s ,M = o, (l) 

where s is the spin of the particular integral. The IM can be given as integrals of the corresponding 
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densities (currents): 

r a+ i + e._id2]. (2) 



j\p. 



Note that for a conformal theory the second term in (J2J is zero. 

An important direction of development is the study of CFT and IQFT with boundary ^3] JH] 
15 . One of the motivations has been the application to string theory and D-branes ^B] |2Q| 
The effect of a particular boundary condition can be given by a boundary state \B), which is a 
vector in the space of states of the bulk theory 

n = I Fj^Fj, 



where Fj and Fj are the representations of the left and right copy of the chiral algebra A labelled 
by j- 

Even if a bulk theory is conformal, it is possible to have a non-conformal, but still integrable 
boundary condition, which permits the survival of infinitely many integrals of motion |13j . Here 
by "survival" we mean that the boundary condition preserves the quantity I s — I s in the sense that 
classically there is no flow of this integral of motion across the edge of the worldsheet. In quantum 
version that means 

(I s -I s )|S) = 0. (3) 

This means that any boundary state can be given as a linear combination of eigenstates \a, j) € Fj 
of the integrals of motion : 

\B) = [ Y,Bi( K )\a,j)®\^jj, (4) 
i a 

where the coefficients in fact depend on the energy scale k, and a labels the eigenstates. 

It is possible to identify between states in Fj <g> Fj and operators in the space of linear maps 
Fj — ► Fj, therefore we can associate a boundary operator B to the boundary state \B): 

i(«)= jY, B iw,j)M- (5) 

J a 
3 

In some cases the boundary operator can be expanded like 

logBO) = logic + V^. (6) 

s 

Thus, if we know a particular family of commuting charges, and find a boundary condition com- 
patible with ©, then we can describe the associated boundary state. Also, it can be inferred that 
the partition function for such a theory, which equals to the overlap of the boundary state with 
the radial quantization vacuum, 

Z{k) = (B(k)\0), (7) 

can be similarly expanded in powers of n, whith coefficients proportional to the vacuum expectation 
values of the commuting IM: 

logZ = logZ + Vii (8) 

s 

where / s = (I s ) arc the vacuum expectation values of the IM and Zq is the partition function of 
the infrared fixed point. 
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Another piece in this mosaic was noticed originally by Dorey and Tateo |y] and corroborated 
later in other works Q [2] 0, namely, that there exists a peculiar relation between integrabe 
field theories and the spectral theory of certain linear ordinary differential equations (ODE). This 
relationship, sometimes called the "ODE/IM correspondence", is still not very clear, but never- 
theless there is substantial evidence for it. In particular, for many integrable field theories it has 
been shown the existence of a special ODE, whose solutions' Wronskian reproduces with amaz- 
ing exactness the functional dependence of overlapping amplitudes of the IQFT on the theory's 
parameters. 

In this context, the problem of study and classification of families of commuting integrals 
of motion {I s } is not uninteresting. As a first step in that direction one can try to find such 
families for the case of free bosonic theory in the bulk. If such ones are found (rather "when" 
than "if" , e.g. we know of the existence of the KdV-Virasoro family of IM, generated by T and 
its composite fields, and others are expected), we can try to identify what possible boundary 
condition is compatible with that family, and find out can an ODE be written in the context of 
the ODE/IM correspondence. Although this task by itself is not so small (unless/until by some 
deeply nontrivial means we gain insight into the problem of general classification of families of IM 
for free bosonic theory), plausible future steps could be the study supersymmetric generalizations 
or massive backgrounds. 

The goal of the present work is to present a method for calculating of commuting integrals of 
motion for free bosonic theory. The only limitation for this method is the problem of solving large 
over-defined square (bilinear) homogeneous algebraic sytems. We present also some results 1 which 
were obtained or verified by applying this method. For N bosonic fields with 0{N) symmetry, this 
includes the (already known) KdV-Virasoro series pQ [22] [2S]> the u O(N)" series recently found in 

and another series, which is related to quantum Boussinesq theory |7j. For the case of two fields 
having Z2 x Z2 symmetry, we obtained a deformation of the O(N) case, which was identified as 
related to the paperclip boundary condition in a previous work [5]. In addition, for the case of iV+1 
fields with symmetry Z2 x O(N), obtained were three new families, which in some sense generalize 
particular cases of the paperclip family. We will also introduce the software implementation of 
the method in question which was developed in the course of this reserch, which could be easily 
adapted to cases of other symmetries as well as other problems involving computation of operator 
product expansions. 

2 Calculating the commuting local IM 

For our purposes it is convenient to consider the theory on a cylinder with complex coordinates 
v, v — x ± iy, where x = x + 2t:R, y € [0, 00). In this case the bosonic integrals of motion will be 
homogeneous. Let us consider the space B of bosonic currents generated by linear combinations 
and normal ordered products (at the same point v) of 



Here and further, for simplicity, we use the notation d = d v . For any current of spin s P s S B we 
can construct the corresponding charge of spin s — 1 : 



C 

where the contour C is taken around the cylinder, and the spin s equals to the number of derivatives 
in P s . We denote the algebra of the charges Q, which is a sum of spin-s subspaces 

It is well-known that free theories have infinite series of commuting integrals of motion which are quadratic in 
the fields, but the aim of this work is the study of other, less trivial families. 



d n X j (v), neN, je {!,... N}. 



(9) 




(10) 
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Figure 1: Semi-infinite cylinder 



Q = 0Q s - 1 . (ii) 

Likewise, B is a sum of subspaces (not subalgebras!): 

oo 

# = 0# s (12) 



8=1 



In fact, the algebra of charges is isomorphic to the space of currents modulo total derivatives: 

(13) 

The same can be stated about the subspaces: 

Q 8 ' 1 = B s /dB s -\ (14) 
The commutator of two charges, I Sl _i € Q Sl_1 and I S2 _i € Q S2_1 , is 

[I Sl _i,I S2 _i] = j> dw j> dzP Sl (z)P S2 (w). (15) 

C w 

It is easy to check that 

[I Sl _i,I S2 _i] = 2m I dw J Sl+S2 _i(w) 
c 

= 27T«i? ;sl+;S2 _2, (16) 



where J Sl+S2 _i(w) e gsi+^-i j s ^h e re sidue of the operator product expansion (OPE) of P Sl £ B Sl 
and P S2 € # S2 , and -R sl + S2 -2 <= Q Sl+;S2 ~ 2 is the corresponding "charge". In this paper I use a not 
very standard normalization of the fields, so the Wick pairing is given by 



dX\z)dX'{w) = - ( rj. (17) 

(z — wy 
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To go to a more physical normalization, like the one in [5], one needs to do the substitution 
— > i\/2X . Likewise, the momenta in the vacuum expectation values in Section 0] and further 



need to be substituted as P 



V2 



P. 



(19) 



It is clear that all charges in Q s 1 can be expressed as linear combinations of basis elements 

/ a b lj = j dzPW ,j = l,..., dimC -1 , (18) 

where P s is an appropriately chosen basis of monomial bosonic currents, which span B s modulo 
total derivatives. For example, let us consider s = 4, and for simplicity assume N = 1 (we have 
only one bosonic field). Then B A can be spanned by 

d 4 x, d 3 xox, (d 2 x 2 ), d 2 x{dx) 2 , {dx)\ 

Taking into account that total derivatives do not contribute to the charges 

d*X = 
d 3 XdX + (d 2 X) 2 = 
3d 2 X(dX) 2 = 0, 

we can perform a Gaussian elimination to obtain a basis modulo total derivatives: 

pj 1] = (d 2 X) 2 . 



(20) 



(21) 



The above is a fairly simple example, but for larger N and s the basis tends to grow rapidly, 
making the problem of finding commuting charges difficult. 



s 


1 


3 


5 


7 


9 


11 


dim(B s+1 ) 


1 


3 


7 


16 


35 


77 


dim(Q s ) 


1 


2 


4 


8 


15 


31 



Table 1. Dimensions of the space of 0(iV)-symmetric bosonic currents for several spins. 

If we want the currents to have some particular symmetry, say O(N) or li2 N , the growth of 
the basis modulo total derivatives will be significantly reduced. 

We can propose the following procedure for calculating of two integrals of motion I Sl and I S2 
of spins si and S2 with symmetry group G(we will denote S3 — si + S2): 

1. For / e {1, 2, 3} , enumerate all possible currents pjf € B Sf , and Pj^+i € B Ss+1 , compatible 
with G, where 3 G {1 . . . m = dimS S/ } and k G {1 . . .n = dimS S/+1 }. Then, by "equating" 
the total derivatives c?Pjy = we obtain a homogeneous linear system of m equations for 
n variables. As we saw in the above example, by Gaussian eliminaion one can choose a 
n — m-dimensional basis(in fact, we must write rank(9/]^ = 0), but it turns out that in the 
cases studied so far the equations were independent) for B Sf+1 . 



2. The integrals of motion are constructed as 



dz 



dim Q B 

E 

h=l 



dim Q s 2 

dw ]T h 2 pl 



(22) 



For each (Zi , Z2 ) G {1 . . . dim Q Sl } x {l...dimQ S2 } calculate the residue of the operator 
product expansion (OPE) of P^+x( z ) P n\i( w )- 
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3. The result of the previous operation will be in B S3+1 . To obtain its value in Q Sa , the 
factorization is done by substituting the solution of the linear system obtained for S3 in step 



4. In order [I S1 ,I S2 ] = 0, the result from step must be zero. By requiring that the coefficient 

before each P^+i to be zero, we obtain a homogeneous bilinear system for and bi 2 . 
Although that system usually seems to be overdetermined(for non-trivial values of s% and 
S2), it turned out that the equations are not independent and solutions were obtained, which 
will be presented in|3| 

Here I will address the question of implementing the above calculations in practice. Like it 
was already mentioned above, the size of the basis in Q s grows rapidly with s, thus making the 
systems of equations for the coefficients of commuting charges prohibitively large not only to solve, 
but even to write down by hand. That's why in order to partially cope with this complexity, it 
is a good idea to develop a computer program, which will handle the ingenuous, though plentiful 
operations related to implementing the consequence of steps described in section^ finally providing 
on its output the bilinear system of equations mentioned afore. Further, one needs to solve that 
essentially quadratic algebraic system. 

For a particular symmetry, the program will have a "combinatorial" module, which enumerates 
all bosonic currents of spin s with that symmetry. We need to invoke the module also for spin 
s + 1. Thus, symbol tables are procuced for spins s and s + 1 of sizes m and n respectively. We 
need to have simple function which calculates the derivative of a spin-s current and expresses it 
as a linear combination of s + 1 currents. Thus we have a linear homogeneous integer-coefficient 
system of m equations for n variables. The solution of the system by Gaussian elimination and the 
further calculations involve rational numbers, that's why in this project I used the GNU Multiple 
Precision Arithmetic Library, gmp. Maybe the most crucial part of the program is the module which 
actually computes the residue of the OPE of given two currents. Essentially, one needs to count the 
contribution of every possible Wick contraction. This is done with standard recursion techniques. 
After that part is built, one needs to put all that inventory together. This software 2 is enough 
flexible, so that by providing a different "combinatorial" module, and with some modifications to 
the OPE module, one might study cases of different symmetry groups or related problems involving 
operator product expansions of currents. 

3 Several Families of Commuting Charges 

3.1 O(N) symmetry 

We consider N free bosonic fields X^, /Lt € {1 . . . N}. The currents can be constructed as polyno- 
mials in terms of the form d ni X ■ d n2 X. It is convenient to use the notation 



1 




(23) 




(24) 



2 The source code is available at http://www.physics.rutgers.edu/~vitchev/ope, 



G 



Correspondingly, possible bases for Q 5 , Q 7 and Q 9 are 



PP = (33) 


P m 


(12)(12) 


pj 1] = (22)(11) 


p[3] 
^6 


= (H)(ll)(ll) 



,[1] 

8 ~ 
,[2] 

8 — 
>[3] _ 



^10 
p[l] 
•VL0 
p[2] 
VL0 
p[3] 
M.0 

pW 

^10 
p[5] 
M-O 
p[6] 
MXI 
p[7] 



(23)(12) 
(13)(13) 
(33)(H) 
(22)(22) 



P 



- — 

r[5] 

! 8 — 
M 

8 — 

[7] _ 
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(44) 

(22)(11)(11) 
(12)(12)(11) 
(11)(11)(11)(11) 



(23) (23) 
(55) 

(44)(11) 
(33)(H)(H) 

(24) (13) 
(23)(12)(11) 
(13)(13)(11) 

(33)(22) 



"pW" 

p[9] 
M.Q 
p[10] 
"VIQ 

P [n] 

p[12] 

^lO 

p[13] 

p[14] 



(22)(22)(11) 
(23)(14) 
(22)(12)(12) 
(14)(14) 
(22)(11)(11)(11) 
(12)(12)(11)(11) 
(11)(11)(11)(11)(11) 



(25) 



(26) 



(27) 



There are several interesting families of commuting charges. Here we will give the results for 
the several first members of each family: 

(l)KdV-Virasoro This series 3 is related to boundary Sine-Gordon model [22j jZH] Q, and can 
also be written in terms of composite fields built of the energy- momentum tensor- see l|38|l - (|45(l . 



dz - 2P 



[o] 



(28) 



dz —(56 + N)P£ 



[0] 



6R 



[i] 



' ; (20 + N)pM 



(29) 



It = 



dz ~ — (2 + TV) (39 + 2/V)pf 1 + ^(364 + 34/V + 7V 2 )P™ 



(62 + N)pW 



45 



(229 + 34iV + N J )P& 



-12R 



[5] 



s --(26 + /V)P f 



[6] 



P 



[7] 



[3] 



657 



(1464 + 597V + JV 2 )Pj' 



(30) 



(2)"0(iV)" This is the family of commuting integrals of motion associated with spherical branes 
found recently in 0. It cannot be expressed completely in terms of T: 

h = I dz - | (2 + N)Pf + P] 1] (31) 
15 = / dz 150 (4 + iV)(59 + 367V ) P 6 01 " |( 4 + We' 1 - y ( 4 + ^) p i 21 + P i 31 ( 32 ) 

3 The first member in this series as well as in the other series in this paper is Ii = § dzT(z), where T is the 
energy-momentum tensor. 
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dz - -|j-(6 + iV)(199 + 88JV)pj 0] + 4f(6 + AT)(692 + 2397V)P™ 

~(6 + TV) (37 + 12iV)p[ 2] - JL ( 6 + TV) (211 + 97JV)if ] 

2 12 
: (6 + iV) (11009+ 11183iV + 2700Ar 2 )F f | 41 (6 + iV)P^ 5] 



77175 v /v ' 7 

-y(6 + A0Pj 6l +Pi 71 (33) 

(3)" $i, 2 " or A\ ' This family has spins 1, 6n ± 1, and is related to a reduction of the quantum 
Boussinesq theory [7j and boundary Bullough-Dodd model [22] |25) : 

I 5 = y dz i (l07 + 17JV + 15y/{N - 25)(N - 1)) pj 0] - 6P] 11 

-- (-85 + 17N + 15V(JV-25)(iV-l)) pJ 2] + P,! 31 (34) 



117 = 



r/2 

[o] 



— (3621 + 4160iV - 3257V 2 + (609 + 3157V) ^{N - 25) (N - 1)) p| 
4 ^ (-3171 - 2810iV + 325iV 2 + (441 + 3157V) yJ(N - 25) (N - 1)) pj 1] 

(l31 + 9iV + 7y/{N -25)(N - 1)) Pfj 2] 
-— (8931 + 2810iV - 325iV 2 + (3157V - 441) y/{N - 25) (AT - 1)) p] 3] 

" 14400 (~ 16929 + U0N ~ 325iv2 + ( 3157V " 2541) ^(TV- 25) (AT - 1)) P f 
-12Pi 5] 



[4] 



3 



[6] 



-- (3 + 9Ar + 7VW^ 25)(JV- 1)J P t 
+PP (35) 
The densities in this series can be written as 



, 181 - 17 AT - 15 J(N - 25) (AT - 1) ,„ s , , , 

P 6 = 8P 3 + ^ ^ '-{dTf (36) 



P 8 = 16P 4 + I ^61 — 9JV — 7-v/(AT- 25)(AT- 1)) 2T(dT) 2 



325N 2 ~ 6590N + 12849 + (315N - 2499) V /(A^ - 25)(A^ - 1), 



480 



T<9 4 T, (37) 



where we use the regular part of the composite fields built from derivatives of the energy- momentum 
tensor: (cf. (JUJ)) 

T=\{\1) (38) 
<9T=(12) (39) 



2T 2 = (11)(11) - 2(22) 



8T 3 = (11)(11)(11) - 6(22)(11) - 12(12)(12) + 3(33) 



(5T) 2 = (12)(12)--(33) 



(40) 
(41) 
(42) 



16T 4 = (11) 4 - 12(22)(11)(11)- 48(12)(12)(11) + 12(33)(11) 
+20(13)(13) - 8(22)(22) - ^(44) 



(43) 



2T(dTf = (12)(12)(11) + 3(23)(12) - |(13)(13) + |(22)(22) - ^(33)(11) + ^(44) (44) 



2Td 4 T = 2(13)(13) - ^(44) - 8(23)(12) - 2(22)(22) 



(45) 



3.2 Z 2 x Z 2 symmetry 

Here we will use the notation 



(m,n 2 , . . . , n a \m lt ma, ...,mp)= d ni Xd n2 X . . . d n «Xd" ll Yd m2 Y . . . d m ?Y. 



(46) 



Where no ambiguity is possible, we will omit the commas. Possible choice of bases for two bosonic 
fields with Z 2 x Z 2 symmetry is 



(47) 





6 
[7] 

[8] 

[9] 

,[10] 



(1122|) 

(1111111) 
(1111111) 
(1111111) 
(1111111) 



(48) 



9 



P m 

pP] 
p[3] 
p[4] 



P 



[■"5] 

3 

p[6] 

pin 

p[8] 
^8 
p[9] 
^8 

[10] 



p, 



[11] 

[12] 



(33|11) 

(44|) 

(13|13) 

(1133|) 

(1223|) 

(144) 
(14|12) 

(|111122) 

(|1133) 

(11|1122) 

(12|1112) 

(H133) 

(22|1111) 



TT3] — — 

i — 
[14] 

i ~ 

[15] = 

[16] = 

>[17] _ 

8 ~~ 

[18] = 

[19] = 

[20] = 

>P1] _ 

8 — 

,[22] _ 

8 — 

>P3] _ 

8 — 

>P4] _ 



(22j22) 

(1111|22) 

(1112|12) 

(23|12) 

(1122|11) 

(|1223) 

(111122|) 

(111111111) 

(11|111111) 

(lllljllll) 

(111111|11) 

(11111111|) 



(49) 



It turns out that the "O(iV)" (when N = 2) family has a Z2 x Z2-symmetric deformation, related 
to the "paperclip" model (SJ. Here are the densities of some of the commuting integrals of motion 
of that family: 



An 2 + 7n + 2 [0] _ 4n 2 + 9n + 4 m 
3(3n + 2)(3n + 4) 6 ~ 3(3n + 2)(3n + 4) 6 

n + 2 pP] , + 2 ) p[3] , g pW 

6(3n + 4) 6 (3n + 2)(3n + 4) 6 6(3n + 2) 6 

(50) 
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96 + 500n + 860n 2 + 575n 3 + 131n 4 
120n(2 + n) ' 
2 + lln + 6n 2 D[ i] 



,[o] 



+ 4 
+ {4+Un + 7n 2 )PP 
32 + 232n + 554n 2 + 473n 3 + 131n 4 [3] 

120n(2 + n) 6 
4+13n + 6n 2 [4] 



4 

(2 + 5n)(4 + 9n + 4n 2 ) [5] 



An ~ 6 



(8 + 5n)(2 + 7n + 4n 2 ) [6] 

4(2 + n) 6 
(2 + n)(2 + 5n)(4 + 5n) [7] 

120n 6 
(2 + n)(2 + 5n) [8] 



n(8 + 5n) [9] 

8 ^ 
(6 + 5n)(8 + 5n) p[10] 

120(2 + n) 



Pf' 1 (51) 



3.3 The Z 2 x O(N) ("cylindrical") case 

Let us consider N + 1 bosonic fields (X, Yi , . . . Yjy) = (X, Y) . There are three non- trivial solutions 
for I 3 and I 5 with Z 2 x O(N), and it seems plausible that these are the first members of infinite 
families of commuting charges. The results can be expressed in terms of only X and the regular 
part of the composite fields of the energy-momentum tensor of Y. Namely, these composite fields 
are 

T Y = -dY ■ dY 
2 

dT Y = dY ■ d 2 Y 

2T Y = {dY ■ dY) 2 - 2d 2 Y • d 2 Y 

8T Y = {dY ■ dY) 3 - 6{d 2 Y ■ d 2 Y){dY ■ dY) - \2(dY ■ d 2 Y) 2 + Zd 3 Y ■ d 3 Y 
{dT) 2 = {dY ■ d 2 ?) 2 - -^d 3 Y ■ d 3 Y . (52) 

In terms of those composite fields, the currents can be written as: 
Solution 1: One solution (Solution 1) is 

P 4 = AT Y + 1 (-9 -3N- y/(N-25)(N -if) {dX) 2 

+ \2{dX) 2 T Y + i (2 + N + y/(N - 25){N - 1)) (dX) 4 (53) 
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(307 + 106 A + 19A 2 + (53 + 13N)J(N-25){N-1)) 
+ ~ 1620 ^ 



4 (l6 + 2A + y/(N -25)(N -1) 
9 

^33 + 3A + y/{N -25)(N 
9 

(l31 - 317V + 8A 2 + (31 + 8N)y/(N-25)(N-l) 



dXd z XdT Y 
(d 2 X) 2 T Y 

N~T)) 

(dx) 2 (d 2 x) 2 



54 

2 f 8 + N + y/(N-25)(N-l)) 

-± „ L (OXfTy 

f 58 - 35A + AN 2 + (17 + 4A) J(N - 25)(N - 1)) 

~ HdXf. (54) 

When N = 1, this solution coincides with the paperclip series with n = — 1. 
A second solution (Solution 2) is 

(409 + 467V + N 2 + (-5 + N) J(N - 25) (N - 1)) 
P 4 = 47f-A _ L^ X f 

(ll + JV +v /(JV-25)(JV-l)) 
+ ^ -(dX) T Y 

(ll7 + 2A + A 2 + (15 + AW(A-25)(A-1)) 
+ A 7^ L ( dX )' ( 55 ) 



8 (20 + N - 2J{N - 25) (N - 1)) (ll + N+ J (A - 25) (N - 1)) 

P ' = ~ 15(2 + N) ' T y + ~ 8 i[aX) ' T 

+ -125 - 7N + 9V(iV - (&Tyf + 4(3x)2r? 



6930 + 1265A + 218A 2 + 11A 3 + (778 + 111 A + 5N 2 )y/(N - 25) (N - 1) 3 ^ 2 
+ 1440(2 + A) 

3 24 v ; y 



^505 + 68 A + 3 A 2 + (27 + 3A) V /(A — 25)(A - 1) 



(9X) 2 (a 2 A) 2 



96 

69 (95 + A 2 + (13 + A) V /(A-25)(A-1)) 



(<9X) 6 . (56) 



480 

For A = 1, this solution coincides with the paperclip solution, but with n = 2. 
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There is also a third solution (Solution 3), namely 

fl3 - 277V + 2N 2 + (5 - 2N)y/(N- 25) {N - l)) 

Pi = ATl + ± J -(d 2 X) 2 

6 

+2 (7 - N + y/(N- 25)(JV- 1)) {dX) 2 T Y + {dX) A (57) 
8 f-5 + 2N+ y/{N -25){N -1)) 



15(2 + TV) 

" ^ " ^(aTy) 2 + 4(dX) 2 T y 



-15 + 3N-y/(N-25){N-l) ffyr , a , 



15 



190 - 61 AT - 7iV 2 + 227V 3 + (_38 _ 417V - UN 2 ) y/ (N - 25) (N - 1) 

180(2 + N) 

2 (9 + 37V - ^(7V-25)(7V-1)) 



d 3 X 2 



7 + 57V- y/(N- 25) (AT - 1) , 2 



(9 2 X) 2 Ty 

1 \ax) 2 (a 2 x) 2 + 2(ax) 4 r y 



40 - 5A^ + A^ 2 + (-8 - A0^(N - 25)(7V - 1) , . s 



6(2 + TV) 

1 (9X) 6 (58) 



-5 + 2AT + ^/(AT - 25) (AT - 1) 



15(2 + AT) 

Like Solution 1, for Af = 1 this one coincides with the "paperclip" solution with n = — 1. 

Summarily, AT = 1, both Solution 1 and Solution 3 coincide with the "paperclip" series with 
n = — 1, and Solution 2 coincides with the "paperclip" series as well, but with n = 2. 

4 Vacuum expectation values 

The vaucuum expectation values of the Z 2 x Z2-symmetric integrals of motion from Section 13.21 
are: 

- 4 + 3n -p4 + 2 + 3n Q4 + p2g2 + 



6(2 + n) 6n 

3 + 2n n2 l + 2n 2 ll + 36n+18n 2 , . 

: P + — a Q + oar, /n r - ( 59 ) 



6(2 + n) 6n 360n(2 + n) 

(6 + 5n)(8 + 5n) 6 _ (2 + 5n) (4 + 5n) 6 
120(2 + n) 2 120n 2 ^ 

8 + 5n_ p4Q2 _2 + 5n p2g2 



8(2 + n) 8n 
(4 + 3n)(8 + 5n) p4 _ (2 + 3n)(2 + 5n) 4 _ 2 + 10n + 5n 2 p2 2 

48(2 + n) 2 48n 2 ^ 8n(2 + n) 

30 + 225n + 250n 2 + 76n 3 2 28 + 137n + 206n 2 + 76n 3 2 
480n(2 + n) 2 480n 2 (2 + n) ^ 

564 + 3410n + 7385n 2 + 5680n 3 + 1420?i 4 
60480n 2 (2 + n) 2 



(60) 
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= n(8 + 7n)(10 + 7n) p8 (2 + n)(4 + 7n)(6 + 7n) 
7 42(2 + n) 2 (2 + 7n) 42n 2 (12 + 7n) W 

, 2,(10 + 7n) , p6Q2 + 2(2 + n)(4 + 7n) p2Q6 + p4Q4 



3(2 + n)(2 + 7n) 3n(12 + 7n) 

7t(5 + 4n)(10 + 7n) i36 (2 + n)(3 + 4n)(4 + 7n) 6 

" 9(2 + n) 2 (2 + 7n) + 9n 2 (12 + 7n) Q 

6 + 49n + 28n 2 . 2 20 + 63n + 28n 2 2 _j 
-P 4 Q 2 H — — — P Q 



3(2 + n)(2 + 7n) 3n(12 + 7n) 

122 + 1407n + 1792n 2 + 612n 3 4 420 + 1583n + 1880?i 2 + 612?i 3 4 
+ 180(2 + n) 2 (2 + 7n) + 180n 2 (12 + In) ^ 

7(40 + 338n + 985n 2 + 816n 3 + 204n 3 ) a 

H i - -; r -P 2 Q 2 



30n(2 + n)(2 + 7n)(12 + 7n) 

1 9,n 3 + 74fi7fm 4 4- 1 , 



2520 + 26166n + 102459n 2 + 136612n 3 + 74676n 4 + 14552n 5 2 



1260n(2 + n) 2 (2 + 7n)(12 + In) 
3720 + 30202n + 89149n 2 + 121284n 3 + 70844n 4 + 14552n 5 



+ 1260n 2 (2 + n)(2 + 7n)(12 + 7?i) 

68760 + 63214 

+ 



68760 + 632142n + 2264647n 2 + 4095840n 3 + 3708040n 4 + 1610448n 5 + 268408n 6 



151200n 2 (2 + n) 2 (2 + 7ra)(12 + In) 

(61) 

5 The vacuum expectation values and the quasiclassical ap- 
proximation 

In this section we will demonstrate the ODE/IM correspondence for the "paperclip" series from 
Section 13.21 The corresponding equation is: 

nP 2 p z f n 4- 9 1 \ 

-d 2 *(z) + k 2 (1 + e?) n 9(z) + — y^VW + [-^-Q 2 + 4 J = °- (62) 

One can try to calculate the Wronskian of l|62|l in the "quasiclassical" limit k — > oo. To this 
end, we use the ansatz 

¥(*) = exp[ K S*_i(z) + S (z) + -Sx(z) + \s 2 (z) + . . .] (63) 

It is possible to show that the expansion of the Wronskian is 

W[f+, /_] = exp[ K< S_i(-oo) + S (z) + -S^-oo) + ^S 3 (-oo) + -W-oo) + ...}. (64) 

K K K 



Calcultions yield: 



, , P 2 Q 2 1 

si-' - t + T + S < 65 > 



(-oo) = n 4 2 + n 4 n(2 + n) , 2 

3 24(2 + 3n) 24(4 + 3n) V 4(2 + 3n) ^ 

n(3 + 2n) d2 (2 + rz) (1 + 2rc) 2 2 + 6ri + 3n 



2 



p 2 v ' ,y ' ' 11 (661 

24(2 + 3n)(4 + 3n) 24(2 + 3n)(4 + 3n) ^ 192(2 + 3n)(4 + 3n) v ; 
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' 40(2 + 5n)(4 + 5n) 40(6 + 5n)(8 + 5n) ^ 

3n 2 (2 + n) 3n(2 + nf 

8(2 + 5n)(4 + 5n)(6 + 5n) W 8(4 + 5n)(6 + 5n)(8 + 5n) ^ 

» 2 (4 + 3n) p4 (2 + n) 2 (2 + 3n) 4 

16(2 + 5n)(4 + 5n)(6 + 5n) 16(4 + 5n)(6 + 5n)(8 + 5n)^ 

3n(2 + n)(2 + 10n + 5n 2 ) 2 2 

8(2 + 5n) (4 + 5n) (6 + 5n) (8 + 5n) ^ 



+ 



n(30 + 225n + 250n 2 + 76n 3 ) p2 
160(2 + 5n)(4 + 5n)(6 + 5n)(8 + 5n) 
(2 + n)(28 + 137n + 206n 2 + 76n 3 ) 2 
+ 160(2 + 5n)(4 + 5n)(6 + 5n)(8 + 5n) ^ 



+ 



24 + 140n + 290n 2 + 220n 3 + 55n 4 
640(2 + 5n)(4 + 5n)(6 + 5n)(8 + 5n) 

15n 3 8 15(2 + nf 



448(2 + 7n)(4 + 7n)(6 + In) 448(8 + 7n)(10 + 7n)(12 + In) W 

15n 3 (2 + n) 

16(2 + 7n)(4 + 7n)(6 + 7n)(8 + 7n) V 

15n(2 + n) 2 - 6 

16(6 + 7n)(8 + 7n)(10 + 7n)(12 + 7n) V 

45n 2 (2 + n) 2 pi . 

32(4 + 7n)(6 + 7n)(8 + 7n)(10 + 7n) W 

5n 3 (5 + 4n) p6 

32(2 + 7n)(4 + 7n)(6 + 7n)(8 + In) 

5(2 + n) 3 (3 + 4n) 6 

32(6 + 7n)(8 + 7n)(10 + 7n)(12 + In) W 

15n 2 (2 + n)(6 + 49n + 28n 2 ) 4 2 

32(2 + 7n)(4 + 7n)(6 + 7n)(8 + 7n)(10 + In) Q 

15n(2 + n) 2 (20 + 63n + 28n 2 ) 2 4 

32(4 + 7n)(6 + 7n)(8 + 7n)(10 + 7ra)(12 + 7n) Q 

n 2 (122 + 1407n+ 1792n 2 + 612n 3 ) pi 
128(2 + 7n)(4 + 7n)(6 + 7n)(8 + 7ra)(10 + In) 

(2 + n) 2 (420 + 1583n + 1880n 2 + 612n 3 ) 4 
128(4 + 7n)(6 + 7n)(8 + 7n)(10 + 7n)(12 + In) Q 

21n(2 + rc)(40 + 338n + 985n 2 + 816n 3 + 204n 4 ) 2 2 
64(2 + 7n)(4 + 7n)(6 + 7n)(8 + 7ra)(10 + 7ra)(12 + In) Q 
n(2520 + 26166n + 102459n 2 + 136612n 3 + 74676n 4 + 14552n 5 ) p2 

896(2 + 7n)(4 + 7n)(6 + 7n)(8 + 7n)(10 + 7n)(12 + 7n) 
(2 + n)(3720 + 30202n + 89149n 2 + 121284n 3 + 70844n 4 + 14552n 5 ) 2 
896(2 + 7n)(4 + 7n)(6 + 7n)(8 + 7ra)(10 + 7rc)(12 + In) Q 
12240 + 110628n + 386834n 2 + 680568n 3 + 606452n 4 + 261786n 5 + 43631n 6 
14336(2 + 7n)(4 + 7n)(6 + 7n)(8 + 7n)(10 + 7n)(12 + In) 
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One can observe that 



Si 
S3 
S 5 
S 7 



where a s are normalization factors, and 

Pi 

03 
05 

turn out to be coefhcients of the expansion 

iogr(2«o = ...+ 



oeih + Pi 
ash + P3 
+ 05 
a 7 I 7 + (3 7 



1 

24 



(69) 



2880 
1 



40320 
1 



215040 



E 

k=l 



B 



2k 



2fc(2fc- l)2 2fe - 1 K 2fc - r 



fhft-l 

where B 2 k are the Bernoulli numbers. 
Thus, for the Wronskian we get 



B2k 



2fc(2fc- 1)2 2 



fc-i • 



W[/ +J /_]~r[2«] exp 



h , h 



(70) 

(71) 
(72) 

(73) 



6 Conclusions 



Commuting integrals of motion are an important feature of integrable field theories, and, in par- 
ticular, conformal field theories which allow integrable boundary interaction. In that case, the 
associated boundary state in such a theory preserves the corresponding family of bulk IM. More- 
over, the boundary state, namely its overlap amplitudes with states of the Fock space of the bulk 
theory, can be expressed in terms of eigenvalues of the integrals of motion [3|[T][2]. In order to 
find families of commuting integrals of motion, associated with such theories, I designed several 
computer programs (available on demand), which are performing the necessary calculation of com- 
mutators of homogeneous bosonic integrals of motion. In this work we showed the first members 
of several families of commuting bosonic integrals of motion, obtained as solutions of the equations 
generated by the programs in question, with symmetries 0(N), Z 2 x Z 2 , and Z 2 x O(N). Some of 
these families are identified as being associated with a particular boundary interaction- the Z 2 x Z 2 
one with the "paperclip" [S] and the O(N) family is related to u O(N) spherical brane model" 0. 
Also, the KdV-Virasoro family is related to boundary Sine-Gordon model (221 QH] Q], and the 
Boussinesq family is related to boundary Bullough-Dodd model [21] For some other families, 
namely the three families with Z 2 x O(N) symmetry, it is not yet very clear what is the boundary 
interaction associated with them. 

Possible future developments might include non-conformal, but still integrable bulk theories, 
like the sausage model |12j . or supersymmetric generalization. 
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